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ABSTRACT 

The possible formation and evolution of dwarf irregular galaxies from material derived from perturbed 
evolved galaxies is addressed via an H I study of a likely example, the peculiar system NGC 5291. This 
system, located in the western outskirts of the cluster Abell 3574, contains the lenticular galaxy NGC 5291 
which is in close proximity to a disturbed companion and is flanked by an extensive complex of numerous 
knots extending roughly 4 f north and 4 f south of the galaxy. In an initial optical and radio study, Longmore 
et al. (1979, MNRAS, 188, 285) showed that these knots have the spectra of vigorous star-forming regions, 
and suggested that some may in fact be young dwarf irregular galaxies. High resolution 21 -cm line 
observations taken with the VLA are presented here and reveal that the H I distribution associated with this 
system encompasses not only the entire N-S complex of optical knots, but also forms an incomplete ring or 
tail that extends approximately 3 ' to the west. The H l associated with NGC 5291 itself shows a high 
velocity range; the Seashell is not detected. The formation mechanism for this unusual system is unclear and 
two models — a large, low-luminosity ram-swept disk, and a ram-swept interaction — are discussed. The H i 
in the system contains numerous concentrations, mostly along the N-S arc of the star-forming complexes, 
which generally coincide w ith one or more optical knots; the larger H I features contain several x 1 0 Mq 
of gas. Each of the knots is compared to a set of criteria designed to determine if these objects are bound 
against their own internal kinetic energy and are tidally stable relative to the host galaxy. An analysis of the 
properties of the H I concentrations surrounding the optical star-forming complexes indicates that at least 
the largest of these is a bound system; it also possesses a stellar component. It is suggested that this object 
is a genuinely young dwarf irregular galaxy that has evolved from the material associated with the system 
and that this entire complex contains several proto- or young dwarf irregular galaxies in various stages of 
development. We are therefore witnessing the early evolution of a number of genuinely young galaxies. 

Given the evident importance of the NGC 5291 system as a “nursery” for young galaxies, careful modeling 
is required if we are to understand this remarkable galaxy. © 1997 American Astronomical Society. 
[S0004-6256(97)02 110-9] 


1. INTRODUCTION 

NGC 5291 lies at a = 13 h 47 m 24?3, 8= - 30° 24' 26" 
(J2000), in the western parts of the cluster Abell 3574 2 
(Abell 1958) of which it is probably a member. We adopt a 
value for H 0 of 75 km s” 1 Mpc -1 ; the duster recession 


'Royal Observatory, Blackford Hill, Edinburgh EH9 3HJ, United Kingdom. 
2 A.k.a. Shapley 1346-30 (Shapley 1936), 1C 4329 (de Vaucouleurs & de 
Vaucouleurs 1964), and Klemola 29 {Klemola 1969). 


velocity of 4350 km s _l (e.g., Richter 1984) then implies a 
distance of 58 Mpc. Pedersen et aL (1978) presented and 
briefly discussed deep images of the system taken with the 
ESO 3.6 m telescope, drawing attention to the immense com- 
plexes of optical knots which extend several galactic diam- 
eters to the north and south and hypothesizing that these are 
H II regions (or at least, sites of recent star formation) from 
their blue colors. In the Second Reference Catalog of Bright 
Galaxies (de Vaucouleurs et al. 1976) NGC 5291 is classi- 
fied as a peculiar elliptical (Ep?) with a nearby companion. 
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Table 1. Properties of NGC 5291. 



Morphological Type a 

SA.0 + 

Magnitude m h b 

15.13 

Optical Dimensions 6 

i: ixo.7 

Adopted Distance c 

58 Mpc 

H i Mass jW H , ad 

4.9x l0'°/W o 


Longmore et ai 1979. 
b NASA IPAC Extragalactic Database. 
c // 0 = 75 km s 1 Mpc" 1 . 
d Single-dish H i mass for entire complex. 

Longmore et ai (1979, hereafter L79), using deep IllaJ 
plates from the UK Schmidt Telescope for the Southern Sky 
Survey, reclassified NGC 5291 as SA.0+. 

Table 1 summarizes the properties of the NGC 5291 sys- 
tem, while Fig. 1 (Plate 74) shows a 5-band image of NGC 
5291 taken by Pedersen et ai (1978) with the ESO 3.6 m 
telescope. The main bodies of the galaxies NGC 5291 and its 
companion (dubbed ‘‘the Seashell” by L79 because of its 
appearance) are themselves relatively unremarkable. The 
overall lenticular morphology of the former is evident in the 
bulge-disc structure, and its late lenticular sub-type in the 
weak, smooth, possibly spiral, structure and dust lane evident 
to the NE. The Seashell is basically smooth in texture, and 
therefore probably also lenticular in underlying type; perhaps 
slightly earlier than NGC 5291, since no dust structure is 
visible. It displays dramatic evidence of its tidal distortion in 
its contorted disk morphology and sharply defined edge. 
Nevertheless, the most unusual feature of this double system 
is not so much the morphology of the main bodies of the 
galaxies but the clearly-associated swarms of knotty conden- 
sations, ranging in size down to the resolution limit of the 
image, and extending ~4' both north and south of the pair. 

Longmore et ai (L79) obtained optical spectra of NGC 
5291, the Seashell and of a number of the outlying knots, 
confirming that these last are indeed the locations of vigor- 
ous star formation as suggested by Pedersen et ai (1978), 
and have the spectra of H ll regions. Longmore et ai also 
made 21 -cm observations of the system, using the 64 m 
Parkes (HPBW=15') and 100 m Effelsberg (HPBW-9') 
telescopes, which revealed that ~5X10 i0 M© (7/ 0 =75 
km s" 1 Mpc ‘j of H I gas was present in the complex. The 
velocity ranges both of the optical features and of the H I are 
> 650 km s _ 1 and are centered on that of NGC 5291, with 
which they are clearly associated. This H I distribution was 
not resolved, however, although the H I centroid could be 
seen to be offset ~L5 to the west of the optical system 
(primarily, it is now apparent from our new data, as a con- 
sequence of the source geometry). 

In their paper, L79 drew attention to the luminosities and 
sizes of the brighter H ll objects in the complex and particu- 
larly pointed out that they appeared to have all the properties 
of dwarf irregular galaxies, and that if these are bound enti- 
ties, they are eventually likely to be lost from the parent 
galaxy and will assume an independent existence as young 
dwarf systems. Two models describing the formation of 
these objects were then postulated, namely, ram-sweeping 
and local compression of a pre-existing, very extended H I 


disk and consequent collapse of the denser concentrations 
into star-forming cloud complexes, or the formation of these 
complexes from clumps of material tidally removed from the 
main bodies of the galaxies by the effects of the interaction. 

Longmore et ai discussed the overall properties of the 
system and noted the very evident interaction between NGC 
5291 and the Seashell. However, they pointed out that the 
high relative velocity of the two galaxies ( — 600 km s~\ in 
an apparently retrograde sense) which they measured made 
the creation of such massive and extensive tidal features 
seem unlikely. The asymmetric, curved optical structure and 
the geometrical offset of the H I centroid, together with the 
system’s location in the northwest edge of a cluster, then led 
them to suspect the presence of ram-sweeping. They specu- 
lated that NGC 5291 initially possessed a very extended, 
tenuous disk of H I, and suggested that this disk is now being 
compressed from the east as the galaxy falls into the Abell 
3574 cluster and encounters an intracluster medium. Insta- 
bilities triggered by this ram-sweeping process would then, 
they suggest, have given rise to the observed star-forming 
complexes, while the sweeping-up and compression of the 
H 1 would give rise to the asymmetric distribution they ob- 
served. 

To investigate further these scenarios, and the intriguing 
possibility that the star-forming objects in the extended 
complex are, indeed, dwarf galaxies in the process of forma- 
tion, high resolution H I observations were taken using the 
Very Large Array (VLA) Radio Telescope 3 and are pre- 
sented here. The observations described in this paper of the 
system, the debris field, and its contents show that it is even 
more complex and unusual than previously suspected. 

In Sec. 2 of this paper, we present the observational pro- 
cedures and data reduction; Sec. 3 contains a description of 
the H i distribution and kinematics, and Sec. 4 discusses the 
analysis and interpretation of the results, including discus- 
sions of both the primordial disk and interaction hypotheses. 
Conclusions and a summary are presented in Sec. 5. 

2. OBSERVATIONS AND DATA REDUCTION 

High spatial resolution spectral line observations of the 
21 -cm emission from the region of NGC 5291 were made in 
1989 by STG and R. Byrne with the Very Large Array radio 
telescope (VLA) in its CnB and DnC hybrid configurations, 
using a 64 channel spectrometer with a total bandwidth of 
6.25 MHz. The resulting channel separation of 97 kHz pro- 
vides velocity resolution of 25.6 km s -1 ; the observations 
were made at a central velocity of 4386 km s ~ 1 . Data were 
collected for two polarizations, right and left circular, which 
were then averaged together to enhance the signal-to-noise 
ratio in the final images. 

The hybrid configurations have a long north arm to main- 
tain resolution at extreme southern declinations where pro- 
jection effects for the beam are large. At the declination of 


3 The VLA is operated by The National Radio Astronomy Observatory; the 
National Radio Astronomy Observatory is a facility of the National Science 
Foundation operated under cooperative agreement by Associated Universi- 
ties, Inc. 
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Table 2. Observational parameters. 



CnB Array 

DnC Array 

C + D Array 

Date 

14 May 89 

13 Oct 89 & l Nov 89 


Total No. of visibilities 

84542 

79479 

164021 

Total bandwidth 

6.25 MHz 

6.25 MHz 

6.25 MHz 

Total No. of channels 

63 

63 

63 

Channel separation 

20.7 km s“ 1 

20.7 km s 1 

20.7 km s-' 

Central velocity 8 

4200 km s 1 

4200 km s * 

4200 km s 1 

Beam (FWHM) 

14"x13" 

50"X39" 

26" X 15" 

rms noise in channel maps 

0.60 mJy/beam 

0.71 mJy/beam 

0.56 mJy/beam 

Equivalent brightness temp. 

2.00 K 

0.22 K 

0.87 K 


a Barycentric velocity. 


NGC 5291 the more extended CnB array provides spatial 
resolution of 14"xl3" while the compact DnC array, with 
lower spatial resolution (50 w X39 r/ ) is more sensitive to weak 
extended emission. A third set of data was produced by com- 
bining that from both arrays — this takes advantage both of 
the higher resolution of the CnB array data and of the higher 
sensitivity of the DnC array observations. The spatial reso- 
lution of the combined data set (hereafter referred to as the 
combined array data, or the C+D array data) is 26"X15". 
The observational parameters for all three sets of data are 
shown in Table 2. 

Calibration and editing of the data were performed using 
the AIPS data reduction package available from NRAO. The 
phase and amplitude response of the receiving system were 
calibrated using the continuum sources 1354 152 and 1313- 
333 for the CnB array observations, and 1313-333 for the 
DnC array observations. The calibrators were observed every 
thirty minutes for the duration of the observations. These 
calibration sources and the instrumental bandpass were in 
turn calibrated for radio frequency flux using the primary 
standard 3C 286. 

The data from each configuration was edited and cali- 
brated independently and converted to separate high resolu- 
tion and high-sensitivity image cubes; the two data sets were 
also combined in the uv -plane to produce the C + D com- 
bined array dataset. 

The data reduction was complicated by the presence of an 
internal interference spike at 1400 MHz, uncomfortably 
close to the observed frequency (1400.8 MHz) of the Hi 
hyperfine transition at the systemic velocity of NGC 5291. 
Calibration involves the use of a continuum data set which is 
normally created by averaging together the central two-thirds 
of the channels in the bandpass. The continuum data set cre- 
ated during our observations was therefore affected by the 
interference spike, which was present in a few of the central 
channels. Rather than attempting to identify and remove the 
contaminated data piecemeal, a new interference-free con- 
tinuum data set was created by identifying 21 interference- 
free channels (those well away from the problematic fre- 
quency) in the line dataset and averaging them in the 
uv -plane. 

After calibration, the line data were examined and edited. 
The channels affected by the interference spike were easily 
identified and the affected data were removed by “clipping'’ 
data points above a specified value determined from an ex- 
amination of the unaffected channels. This process mini- 


mized the removal of uncontaminated data points and hence 
resulted in the removal of very few data points ( — 2% of the 
total) in unaffected channels. The contaminated channels, 
however, suffered the loss of approximately half of the data, 
thus degrading their signal-to-noise. Fortunately, not more 
than three channels were severely affected, due to the narrow 
frequency range of the interference spike. The poorer quality 
of the data from the heavily clipped channels resulted in 
some loss of signal during the blanking and integrating pro- 
cesses (described below), and slightly reduced the column 
densities and H I masses calculated from the integrated 
dataset. 

The uv data set was Fourier transformed to produce im- 
ages of the intensity distribution of the source for each chan- 
nel. The data cube of channel maps resulting from this pro- 
cess includes not only the 21 -cm line emission from the 
observed source, but also the continuum emission from other 
sources. Continuum estimation was performed by combining 
channels which were free of line emission from both ends of 
the data cube, thirteen from the low frequency end and 
eleven from the high frequency end. This averaged con- 
tinuum map was then subtracted from each channel of the 
cube. 

Since the images are produced using a fast Fourier trans- 
form, any unsampled points in the uv -plane are automati- 
cally set to zero. In addition, the finite nature of the uv data 
set produces a truncation effect, so that the fast Fourier trans- 
form process results in a synthesized beam with extensive 
sidelobes: the “dirty beam.” The continuum-subtracted data 
cube, consisting only of 21 -cm line emission from the 
source, is actually a representation of the source brightness 
distribution convolved with the dirty beam. Confusion 
caused by the dirty beam was corrected by using the AIPS 
CLEAN routine (Hogbom 1974; Clark 1980). The resulting 
point-source model was then convolved with the “clean 
beam,” a 2D Gaussian with the same full-width half- 
maximum (FWHM) as the dirty beam, to remove grating 
responses and sidelobes. 

The data cubes were CLEANed down to 2X the average 
rms noise (<r) as determined by a statistical analysis of the 
signal-free channel maps. Synthesized beamsizes represent- 
ing the effective spatial resolution of the Gaussian CLEAN 
beams were calculated for all three data sets; these values are 
provided in Table 2. The average rms noise after CLEANing 
was determined to be 0.60 mJy per beam solid angle for the 
CnB array data, 0.71 mJy per beam solid angle for the DnC 
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RIGHT ASCENSION (B1950) 


Fig. 2. Grayscale image of H I integrated intensity from the CnB array data, with 2<r (0.36X 10 20 atom cm -2 ) contour from DnC array data. CnB and DnC 
array beamsizes indicated by lower (14"X 13") and upper (50"X39") ellipses in lower left comer. Crosses indicate the optical knots identified by L79 (A - H) 
and the H [ features discovered with the VLA (l-M). 


array data, and 0.56 for the combined array data (Table 2). 

After the data cubes were CLEANed, final imaging was 
performed by the integration of signal from —40 contribut- 
ing channels. This was done using the AIPS routine momnt, 
which conditionally blanks the data cube using a mask cre- 
ated by smoothing the data in both velocity and space (cx,S) 
and then integrates the blanked dataset. Blanking is per- 
formed in an attempt to isolate the signal from the noise by 
suppressing the noise. The data for NGC 5291 were blanked 
and integrated using a mask created by Hanning smoothing 
in velocity, and convolving in right ascension and declina- 
tion with a Gaussian corresponding to twice the FWHM of 
the clean beam. The blanking and integrating levels were set 
to approximately 1.5 times the rms noise level on the un- 
smoothed data. Moment maps representing the H I column 
density (integrated flux) distribution (0th moment), velocity 
field (1st moment), and velocity dispersion (2nd moment) 
were created by this process. 

3. RESULTS 

3.1 H l Spatial Distribution 

Moment maps from the high sensitivity DnC array data, 
the high resolution CnB array data, and the combined CH-D 
data were produced using the method described above. Fig- 
ure 2 shows a grayscale image of the integrated H 1 flux 
distribution from the CnB array data overlaid with the 20- 
contour level of the H I distribution from the DnC array data. 


We have indicated the positions of the original optical knots 
identified by L79 (A-H), as well as the discrete H I features 
(I-M) detected in the VLA observations that also appear to 
be associated with optical knots identified from either the 
optical image in Fig. 1 or the Digitized Sky Survey image. 3 4 
Also identified are the locations of NGC 5291 and of the 
Seashell. The H / forms what appears to be an incomplete 
ring structure extending about 9\9 north-south (163 kpc) and 
5.6 east-west (93 kpc). 

The single-dish observations by L79 found the centroid of 
the H I to be offset — L5 to the west of the optical objects. 
This appears to be an effect of the presence of the westward 
arc, which was unresolved by the 9 r and 15' beams of the 
Effelsberg and Parkes telescopes. To test this, we smoothed 
the integrated intensity image from the present DnC array 
observations with a 9' Gaussian to replicate the resolution of 
the Effelsberg observations. This indeed produced a smooth 
distribution with the peak intensity located — V.5 to the west 
of NGC 5291. 

A grayscale image of the H I integrated flux from the 
combined C + D array data is shown in Fig. 3. The optical 
image from the Digitized Sky Survey (derived from a UK 
Schmidt IllaJ plate) overlaid with these H I contours is 
shown in Fig. 4. The H I is highly structured on scales at 
least down to the CnB beam size and many of the concen- 


4 Obtained using Sky view; see acknowledgments. 
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Fig. 3. Grayscale image of the H i integrated intensity from the C + D combined array data. Contour levels are 1.1 (2rr), 5 10, 15, and 20X 10 2n atom cm 
The beamsize (26"X15") is indicated by the ellipse in the lower right corner. Optical and H i knots are indicated with labeled crosses. 


trations are clearly associated with optical knots. 

The spatial association of the H 1 concentrations with the 
larger optical knots is very striking. Figure 3 shows that the 
H I distribution in the complex is generally clumpy; the 
amount and location of H I associated with the optical knots 
observed by L79 are clearly visible in Fig. 4. The greatest 
concentration of H I is associated not with NGC 5291 itself, 
but with optical Knot B. More detailed inter-comparison of 
the optical and H I images shows that many of the H 1 fea- 
tures and clumps have underlying optical knots. Several ot 
these — those associated with the H 1 lying north of Knot A 
and also to the southwest of Knot B, and possibly with Knot 
M as well — have been observed spectroscopically 5 and con- 
firmed to be emission line regions (H II complexes) like 
Knots A to H. The observed H H regions have optical veloci- 
ties in good agreement with those of the associated H l fea- 
tures (Hawarden et ai 1997), while the rest are very similar 
in appearance in Fig. 1 to the confirmed H n complexes. 

The existence of the westernmost H I emission, including 
the concentrations labeled as knots K and L, was not previ- 
ously suspected; examination of the optical image (Fig. 1) 
reveals faint optical knots coinciding with the H J features 
I-M as well. Knots I, J, and M appear to be associated with 
several faint optical features rather than one large one as for 
the other knots. However, we note that these objects could be 
chance superpositions of faint background galaxies; spec- 

5 By the late Dr. David Allen at the AAT. 


troscopy is needed to verify these associations. 

A strong, highly resolved H I feature (clearly visible in 
Fig. 3) lies — 015 west of a point about a third of the way 
from the nucleus of NGC 5291 to the brightest northern knot, 
B. There is no single bright optical feature underlying this 
H 1 concentration, but it does appear to be associated with a 
swarm of faint optical knots in Fig. I, which also have the 
appearance of H II complexes. 

Some H l, - 20% of the total amount in the system, is 
evidently associated with the main body of NGC 5291, most 
of it lying slightly to the west and north of the center of the 
galaxy (Fig. 3). The association of this material with the 
galaxy is confirmed by the correspondence of its flux- 
weighted systemic velocity, 4370 km s _l (see Sec. 3.3, with 
that of the galaxy as determined from optical observations 
(4335 km s _i , L79; 4330 km s _1 , Hawarden et ai 1997). 

While some H l also lies in a similar projected position 
relative to the Seashell, this material has a velocity —4380 
km s~\ several hundred km s _I too high to be associated 
with that galaxy (which has cz=3950 km s -1 : Due 1995 or 
3730 km s^ 1 : L79; see below). The velocity range covered 
in the VLA observations is —3590 to 4870 km s 1 and in- 
cludes both these values; however no H I was detected near 
these velocities. Evidently the Seashell contains little or no 
H l at present; we derive (see below) an upper limit to its H I 
mass (single-channel 5<x detection) of 5.7X 10 7 A#©. 

Perhaps the Seashell was an H l-poor lenticular, a classi- 
fication which is consistent with its relatively smooth (if dis- 
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Fig. 4, Contours of the H I integrated intensity from the C + D combined array data superimposed over the digitized optical image from the UK Schmidt IllaJ 
Sky Survey plate. Contour levels are 1.1 (2a), 5, 10, 15, and 2()X H) 20 atom cm 2 . The beamsize (26"x 15”) is indicated by the ellipse in the lower right comer. 
Optical and H i knots are indicated with labeled crosses. 


torted) structure (even on the UV image of L79) and also its 
moderately red colors (Daly et ai 1987). If any H I it may 
initially have possessed was located in the outskirts of the 
system it would probably have been lost in the encounter 
with NGC 5291, but the apparent absence, anywhere in this 
system, of detectable H l at the velocity of the Seashell 
strongly suggests that this was not the case. Using the cooled 
grating spectrometer CGS4 on the 3.8 m UK Infrared Tele- 
scope (UKIRT), we have secured a /Gband spectrum cover- 
ing the first overtone (A = 2.3 /im) bandhead of CO in stellar 
atmospheres. This spectral feature is unusually weak in the 
Seashell, which implies that it cannot recently have housed a 
nuclear starburst since strong nuclear CO absorption in the 
atmospheres of red supergiants is to be expected in the after- 
math of such an event. Together with its weak optical emis- 
sion lines, this further argues against there having been a 
significant amount of H I in the central parts of the Seashell. 

The galaxy located in the southwest part of the overall H I 
complex is number 435 in Richter’s (1984) list of galaxies in 
the cluster and surrounding area. He classifies it Sc; from 
Fig. 1 we refine this to SA(s)c. Richter gives a heliocentric 
velocity of 4088 km s^ 1 , strongly suggesting that it is a 
member of the cluster. While its H I component is well de- 
tected on our images among the gas associated with the 
southern end of the NGC 5291 complex, its velocity is quite 
different, so that no confusion arises. Evidently Richter 435 
has no connection with the NGC 5291 complex. 

Five arcminutes to the northeast a second prominent gal- 
axy is also included in the VLA primary beam but is not 


detected in H i. This is a barred spiral listed as number 446 
by Richter, who gives a radial velocity of 4409 km s~ 1 , in 
excellent agreement with the L79 velocity of 4420 kms -1 
and a more recent AAT velocity of 4414 km s~ 1 ; its velocity 
is clearly close to that of NGC 5291. Richter classifies this 
system Sd, which suggests that it should be H l-rich and 
should have been detected in our images. However, scrutiny 
of Fig. 1 show's a generally smooth structure, except for 
some curious lateral spurs attached to the bar, suggesting that 
R446 is actually a lenticular or very early spiral: we propose 
a classification R'SB(s)0/a. Such a system would be quite 
likely to contain little H I. 

Figure 5 shows the radio continuum emission as deter- 
mined from those channels free of line emission. There is 
continuum emission associated with the optical body of 
NGC 5291 as well as that of Knot B which is presumably 
due to ionization from O and B stars. The low level (3o") of 
the detected emission precludes quantitative analysis. 

3.2 H l Column Densities and Masses 

Figures 2 and 3 represent the integrated intensity (column 
density) of the H I flux from the system. The column density 
of the H I along the line of sight for a given position (a,S) 
may be determined for a source by integrating the brightness 
over all velocities using the following expression (Mihalas & 
Binney 1981): 

N Hl (ct,8)= 1 . 82 X I0 18 f T H (a,8)dv, 


(1) 
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Fig. 5. Radio continuum contours over optical image. Contour levels correspond lo 3. 5. and Hkr (l<r — 2.7X 10 4 Jy/beam). The beamsize |26 W X15') is 
indicated by the ellipse in the lower right corner. Optical and H I knots are indicated with labeled crosses. 


where 7 B is the brightness temperature in Kelvin and v is the 
velocity in km s” 1 . The peak column densities found for the 
identified features (the galaxies and knots) on the moment 
maps are listed in Table 3. These values were determined 
using the DnC array data, which is most sensitive to ex- 
tended emission. The highest column density in the entire 
complex, 12.6X 1() 20 atoms cm” 2 , is associated with Knot B. 
The peak column densities for NGC 5291 and the H 1 Con- 
centration are 10.7 x 10 20 and 10.6 X 10 20 atoms cm ” 2 , re- 
spectively. Peak column densities for the rest of the labelled 
knots range from 1.4 X 10 20 to 7.4 X 10 2() atoms cm” 2 . The 
detection threshold (2 cr level) is 1 X 10 20 atoms cm -2 . 


The integrated flux from the H I gas may be found by 
integrating Eq. (1) over the solid angle subtended by the 
feature, and expressing the brightness temperature in terms 
of the observed flux S(u). The H I mass may then be calcu- 
lated as follows (Mihalas & Binney 1981): 

A/ Hi =2.36X 10 s D 2 J S(v)dv, (2) 

where , is expressed in solar masses, D is the distance in 
Mpc, S( v ) is the flux in Janskys summed over the area of the 
feature, and v is velocity. The areas used to determine the 


Table 3. Stability analysis of NGC 5291 complex components. 


Object 

(I0 2 ° cm 2 ) 

Ay >0 
(km s ') 

H i radius 
(kpc) 

M h 

(10 ' M ) 

M m i 

( 10 <J A/ 1 

M m , !M h 

Rc, 

(kpc) 

r T 

(kpc) 

NGC 529 1 

10.6 

464.4 

9. 1 J 

773.5 b 

1.60 




Knot B 

12.6 

121.9 

2.3 

3.02 

2.48 

0.82 

61.4 

6.3 

Knot C 

4.4 

306.4 

2.3 

19.08 

0.43 

0.02 

29.6 

1.7 

Knot D 

5.9 

301.2 

1.9 

14.83 

0.54 

0.04 

41.4 

2.5 

Knot F-G 

7.4 

130.5 

1.9 

2.78 

0.91 

0.33 

57.9 

4.3 

Knot H 

5.9 

141.3 

2.3 

4.06 

0.60 

0.15 

74.6 

4.8 

Knot I 

1.4 

87.9 

2.3 

1.57 

0.07 

0.05 

82.1 

2.6 

Knot J 

2.2 

173.6 

1.9 

4.93 

0.23 

0.05 

85.3 

3.9 

Knot K 

3.2 

104.5 

1.9 

1.83 

0.28 

0.15 

81.3 

4.0 

Knot L 

1.6 

71.6 

2.3 

1.04 

0.10 

0.09 

77.7 

2.7 

Knot M 

2.2 

121.9 

2.3 

3.02 

0.13 

0.04 

46.3 

1.8 

H i Concentration 

10.6 

289.2 

1.9 

14.04 

2.03 

0.15 

18.5 

1.8 


“Optical radius (NASA IPAC Extragalactic Database). 

b Keplerian mass M K from Eq. (5); includes a factor of 5 for a dark halo. 
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Fig. 6. Grayscale image of H i integrated intensity from the C + D combined array data. Outer contour indicates 2a level from DnC array data. Inner contours 
indicate boundaries used for determining H i masses of knots. 


masses for each knot are indicated in Fig. 6, and the corre- 
sponding calculated H I masses are listed in Table 3. 

The H 1 masses associated with Knot B and the H I Con- 
centration are roughly comparable: 2.48 X 10 9 A/ Q and 
2.03X 1() 9 A/q , respectively. The Hi mass associated with 
NGC 5291 is 1.60X 10 9 Mq . This value represents only 21% 
of the total H I mass of the complex (see below); over three- 
quarters of the total H l mass is located in the extensions. 
The H I masses of the other knots range from 7 X 10 7 A/ G to 
9X 10 8 A/ G (Table 3). 

The total H I mass detected with the VLA for the entire 
complex is 1.80X10 ]0 A/q . This is approximately 36% of 
that detected in the single-dish observations by L79 (Table 
1). Our estimate does not include the H I mass of Abell 3574 
0435 (9.92X 10 8 A/ o ) which must have been included in the 
large beams of the telescopes at Parkes (15' FWHM) and 
Effelsberg (9 ') used by L79. This would only add 5% to the 
total H I mass, however, and so cannot account for the dis- 
crepancy. 

It is well known that single-dish spectra are prone to sub- 
stantial instrumental effects that make them unreliable for 
weak extended emission-line spectra. The Effelsberg data 
may be an overestimate of the H i present. Employing a very 
wide bandwidth, L79 were observing for broad and weak 
spectral features, which is notoriously difficult with an auto- 
correlation receiver. Indeed the data taken by the Parkes ra- 
dio telescope were not used for this reason. Baseline curva- 
tures can add to or subtract from apparent line strengths, 
especially on very broad lines such as that from the NGC 


5291 complex. It is therefore possible that the discrepancy 
between the masses determined from the VLA and single- 
dish observations is somewhat less than reported. It is ques- 
tionable, however, whether the entire difference between the 
VLA and single-dish fluxes could be attributable to baseline 
errors. Examination of the H I spectra in L79 suggests that an 
overestimate of around 40% is the most that could reason- 
ably be attributed to an unrecognized baseline rise around the 
systemic velocity of 4400 km s -1 . New, high-sensitivity 
single-dish observations, secured with a system having stable 
baselines and a beam at least as large as that set by the VLA 
minimum spacing, would clearly be most valuable to resolve 
this problem. 

There are four possibilities for accounting for the “miss- 
ing H I,” more than one of which could be simultaneously 
true. The single-dish measure could indeed be a 30% to 40% 
overestimate, as discussed above. It is also possible that 
some of the H I could be so extended in velocity and/or space 
as to fall below the surface density detection threshold of the 
observations. We note, however, that if the “missing mass” 
were to fill in the ring, then the H I surface density would be 
approximately 2. OX 10 21 cm" 2 , which is well above the de- 
tection threshold of the DnC array. 

3.3 Kinematics 

Individual velocity channel maps of the D array data in a 
15'X 15' area are presented in Fig. 7. It is clear that H I in the 
NGC 5291 complex has an overall ordered velocity struc- 
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Fig. 7. Individual channel maps from the DnC array data. Contour levels are 2, 5, 10, 15, and 20 a. The heliocentric velocity tor each channel is listed in upper 
right comer. 


ture, with velocities varying (for the most part) in a regular 
manner from low values in the north to high in the south. 

The cumulative velocity field of the gas can be deter- 
mined by calculating the intensity-weighted mean velocity of 
the gas at each position (a,S). This is accomplished by tak- 
ing the first moment of the column density relationship [Eq. 
(1)] with respect to velocity as follows: 


(v(a t S)) 


S)dv 

r-Ma,S)dv 


(3) 


Figure 8 shows the first moment map of th$ complex. The 
isovelocity contours are separated by 25 km s -1 . The wide 
N-S velocity range is evident. The maximum velocity, in the 
south of the complex, is 4750 km s~ 1 ; the minimum, in the 
north, is 4040 km s“ *, giving a (large) total velocity range of 
710 km s" 1 . The flux-weighted heliocentric velocity of the 
gas associated with NGC 5291 itself is 4370 km s - 1 whereas 
that for the entire complex is 4406 km s _l . The difference in 
the flux-weighted velocities show that the gas density is not 
symmetrically distributed in velocity space with respect to 
the galaxy. 

The velocity structure, though necessarily simplified in 
this figure by averaging, nevertheless shows considerable 


complexity just north and south of the main body of NGC 
5291: here the isovelocity contours are much closer-spaced 
than in the outlying regions, and considerable E-W velocity 
shear across, e.g., the H I Concentration of Fig. 2 is also 
apparent. The overall behavior is not very dissimilar to that 
which would be expected from a galaxy seen edge-on, but 
with additional complexities. 

The general velocity trend apparent in the main, eastern 
arc of H I is also visible in the western arc as well, but here 
the gradient appears roughly linear overall. The velocity of 
the (northern) knot L in the western arc is —145 km s -1 
higher than the H I in the eastern arc at the same N-S posi- 
tion. In the southern part of the western arc, an opposite 
trend is (barely) apparent: knots J and K have (average) ve- 
locities — 60 and 50 km s~ 1 lower than the H 1 in the corre- 
sponding parts of the eastern arc. 

The spatial distribution of the velocity dispersion of the 
H I is shown in Fig. 9. The velocity dispersions are calcu- 
lated for each position by taking the second moment of Eq. 
(1) with respect to velocity: 


\(v 2 (a,S)) = 


lr_ y: T H (a,S)(v(a,S)-(v)) 2 dv 
V f%T B (a,S)dv 


(4) 
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It can be seen in Fig. 9 that there are two regions of very 
high dispersion, on the order of 100 km s -1 . These regions, 
one encompassing both NGC 5291 and the H I Concentra- 
tion, and one associated with Knot D, imply the presence in 
the line of sight of bodies of gas with widely ranging veloci- 
ties (see Fig. 1 1). These are further discussed below. 

4. DISCUSSION 

4.1 Physical Nature of the System 

The obvious, strong, distortion of the Seashell, presum- 
ably by interaction with NGC 5291, immediately suggests 
that the far-flung outlying structures in the system were 
caused by tidal effects from this interaction. However, the 
large velocity difference between the two galaxies measured 
by L79 (640 km s _l ), and its apparently retrograde sense, 
made the great amplitude of the resulting disruption difficult 
to understand in terms of tidal effects. This led them to seek 
alternative hypotheses for the origin of the outlying features. 

Since the major problem with the tidal disruption scenario 
was the difficulty of moving so large a mass of material to 
such large radii, they hypothesized instead a pre-existing, 
very extensive, and presumably primordial H I disk. This 
was a logical extrapolation of the properties of several other 


H i-rich early-type systems listed by L79, and analogous to 
the then-undiscovered Malin 1 (Bothun et al. 1987; see also 
Sprayberry et al. 1995 for more details of the properties of 
this class of galaxies). 

Longmore et al. then proposed two trigger mechanisms 
for the intense but localized star formation activity seen in 
the outlying complexes. The first suggestion was that NGC 
5291 has suffered an almost orthogonal, off-center, penetra- 
tion by a companion galaxy (though probably not the Sea- 
shell) which has caused an eccentric kinematic density wave 
to propagate through the disk, producing a ring galaxy with 
offset nucleus (see, e.g., Lynds & Toomre 1976; Theys & 
Spiegel 1976). They tentatively identified Richter 446 (to the 
northeast) as a potential culprit, but did not pursue the hy- 
pothesis because the ring was incomplete. 

The second trigger mechanism was suggested to L79 by 
the location of NGC 5291 in the outskirts of a galaxy cluster, 
the concentration of the H II complexes in an arc convex 
toward the cluster and centered on the galaxy and by the 
westward offset of the centroid of the H I emission. These 
led L79 to speculate that the postulated H I disk is being 
ram-swept, more or less from one edge, as the galaxy moves 
mostly eastwards through an intracluster medium, and that it 
is this process which has triggered star formation, by com- 




1437 


MALPHRUS ET AL: NGC 5291 


1437 


4000 4200 4400 4600 4800 



13 44 45 30 15 00 

RIGHT ASCENSION (B1950) 


Fig. 8. Velocity (first moment) map from the C + D combined array data. Isovelocity contours are 25 km s 1 apart and run from 40(X) km s to 4725 km s 


pression of the previously too-tenuous H I disk. 

We have noted that the H I offset, at least, now appears to 
reflect the detailed distribution of the gas as seen at low 
resolution. The kinematics, but especially the morphology of 
the H I in this system as revealed by the higher resolution 
VLA data presented here, encourage us to reinstate as a pos- 
sible scenario the original interpretation of the outlying fea- 
tures as the products of tidal disruption consequent on an 
encounter, most likely with the Seashell. The ring-arc seen in 
the H 1 may represent tidally removed debris along the path 
of motion of the component galaxies of the system, similar to 
the Magellanic stream. Alternatively, the arms, counter-arms 
and bridges of an interaction, first discussed by Toomre & 
Toomre (1972), are seen at very large projection. 

We also note that the western H 1 arc repairs the “incom- 
plete ring” noted by L79 as evidence against a possible or- 
thogonal penetration by some galaxy other than the Seashell. 
As they suggested, Richter 446 is the most probable candi- 
date as it does appear poor in gas and possibly disturbed, as 
might befit a system which has passed though, and been 
swept clean by, a giant neighbor. 

In the rest of this section we discuss the case for ram- 
sweeping, and revisit the “primordial disk” and “tidal dis- 
ruption” scenarios for the creation of the outlying complexes 
in the NGC 5291 system, and thus of the candidate “young 
galaxies” lying within them. There is supporting evidence 
for both creation models, as well as for ram-sweeping. Each 
of the two formation scenarios, however, has limitations 


based on analysis of the existing optical, radio, and x-ray 
data. 

4.1. 1 Ram sweeping 

As NGC 5291 is located in the northwest part of the Abell 
3574 cluster, ram-sweeping would be expected to produce 
compression along the leading (east) edge, assuming that the 
system is moving toward the center of the cluster, located to 
the southeast. On all hypotheses one would expect to find 
that the star formation regions (SFRs) are associated with H 1 
concentrations. If ram-sweeping is occurring, one further ex- 
pects that since the SFRs are presumably of higher density 
than the less concentrated H I, the latter will tend to be dis- 
placed in tail-like structures trailing “downwind” from the 
SFRs. Furthermore, one expects the density gradient of the 
H I to be steeper on the side towards the compressing wind 
than on the opposite side. 

The gradient in the H 1 distribution is indeed higher on the 
eastward edge of the east side of the ring-arc (Fig. 4), which 
provides some evidence that compression of the eastward 
edge of the H I is occurring, at least on the east side of the 
complex. The distribution of the H I around the H II com- 
plexes (SFRs) B and H, in particular, exactly fulfills the pre- 
dictions of this scenario for ram-sweeping by a flow from the 
east: the densest H I lies on or slightly west of the optical 
features, the H I density gradient is steeper to the east and 
shallow, and the H l much more extended, to the west. These 
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Fig. 9. Velocity dispersion (second moment) map from the C + D combined array data. Contours are 10, 25, 30, 50, 75, and 100 km s 


effects are also present in the complex south of knot B: to the 
east a “prow” of H I wraps steep density gradients close 
around a knotty optical feature but tapers more gradually to 
the west. However, we note that this sweeping to the west 
does not appear to be true for the knots in the western arc 
(particularly K and L), which show a slightly higher density 
gradient on their western edges. 

Further evidence that ram-sweeping is operating may be 
indicated by the slightly asymmetric distribution of the 
atomic hydrogen in the vicinity of the optical body of the 
galaxy. As seen in Fig. 4, the gas appears to extend farther to 
the west of the optical emission than to the east, which is 
expected if the H l is being swept to the west. 

To examine further the validity of ram-sweeping, Hawar- 
den et ai 1997 have obtained x-ray images of NGC 5291, 
first with the EINSTEIN Observatory’s Imaging Proportional 
Counter (IPC) in 1980 and then with the tfOSATHigh Reso- 
lution Imager (HRI) in 1996. An archival observation of the 
super-Seyfert galaxy IC 4329 A with the ROSAT Position- 
Sensitive Proportional Counter (PSPC) (now defunct) is also 
available. It was expected that if ram-sweeping was indeed 
occurring to a significant extent, density enhancement (and 
possibly, if the galaxy is moving supersonically, heating) of 
the ICM “upstream” of the Hi should result in enhanced 
x-ray emission from just ahead of the H I being swept. Pos- 
sible ablation effects were also anticipated. 

In fact NGC 5291 is a quite prominent extended x-ray 
source. Although the resolution (2!6) is not sufficient to 
examine the distribution of the x-ray emission in detail, the 


EINSTEIN image (Fig. 10) shows a source coincident with 
the galaxy and possibly extended N-S along the H II com- 
plexes and an apparent “tail” to the west, as might be asso- 
ciated with ablation by a westward flow. The north-south 
extent of the x-ray emission is approximately 7', somewhat 
smaller than the 10' extent of the H I gas. Compression of 
the x-ray emission along the eastern edge is also evident. 

The ROSAT images (Hawarden & Chaytor 1997; Ha war- 
den et ai 1997) have confirmed some, but not all, of this 
impression. While there certainly exists extended emission 
which may be spatially associated with the H II complexes to 
the north of NGC 5291, several (variable) sources underly 
the “tail” seen in the EINSTEIN image; when these are re- 
moved the tail is reduced. The ROSAT images lie to the west 
of the H I and H II complexes to the south of the galaxy. 

4.1.2 Large H I disk model 

Longmore et ai concluded from both the optical and H I 
velocity field that the complex of knots and the H I itself 
were entirely associated with NGC 5291. As noted above, 
their interpretation therefore was that of a low luminosity SO 
galaxy surrounded by a large (160 kpc) primordial H I disk, 
similar to the then-unknown Malin 1. Their estimate of the 
total mass in the NGC 5291 system, assuming it was a disk 
in keplerian rotation, was 5:9 X 10 n Afo, making it one of 
the most massive systems then known. 

As set out above, this scenario supposes that the outlying 
star-forming regions in NGC 5291 have been produced by 




1439 MALPHRUS ET AL: NGC 5291 


1439 



Fig. 10. EINSTEIN X-ray image courtesy T. Hawarden. The very strong 
source in the northeast (upper left) of the image is the Seyfert galaxy 1C 
4329A; the NGC 5291 complex is the extended source in the center of the 
frame. North is up; east is left. The grid lines are 2 minutes by 20'. 


local compression of the H I disk by the ram-sweeping pro- 
cess, giving rise to instabilities and collapse of large bodies 
of gas to form the H il complexes now observed in the opti- 
cal. 

The overall kinematics of the gas in the system, as seen in 
Fig. 8, could be consistent with a heavily disturbed massive 
rotating disk centered upon NGC 5291. The eastern gradients 
seen in both the H I and x-ray data, as well as the x-ray 
“tail” to the west are also consistent with a ram-swept disk 
system. 

The previously unsuspected ring-arc morphology of the 
H I complex may argue against the existence of a large disk 
however. As discussed in Sec. 3.2, above, no atomic hydro- 
gen was found “filling in” the ring-arc. The association of 
the westernmost knots ( I— L) is also difficult to understand if 
the NGC 5291 complex is to be interpreted as a large disk. 
Compression and subsequent knot formation along the trail- 
ing edge are then difficult to explain, although it is possible 
that the western arc traces out a (large!) spiral feature. 

There is also evidence that the gas associated with the 
NGC 5291 system is not primordial. Metallicity measure- 
ments of the gas in the knots, approximately one-third solar 
(Due 1995; Hawarden et al 1997), is higher than expected 
for primordial material, and indicates enrichment from prior 
star formation. Either the knots formed from pre-enriched 
material such as found in the outer parts of actively star- 
forming spiral galaxies, or there has been self-enrichment 
occurring in the vigorously star-forming optical knots. This 
last would increase the metallicity of the gas to an extent 
dependent upon the amount of mixing in the object (i.e., 
from stellar winds, etc). 


4.1.3 Interaction model 

Although in many ways the NGC 5291 complex does not 
closely resemble a “classical” interacting system, many 
properties of the extensions could be consistent with the sce- 
nario attributing their properties to a tidal origin. The most 
striking evidence is the optical morphology of the system. As 
stated above, the Seashell, with its strongly distorted disk 
and sharply-defined, curved southwestern edge is evidence 
that an interaction could be of importance, since one is cer- 
tainly affecting (or has recently affected) this galaxy. Its pro- 
jected position close to NGC 5291 makes the latter the logi- 
cal candidate to have caused the disruption, and thereby to 
have been disrupted in its turn to produce the outlying fea- 
tures (though we cannot use these to support the identity of 
NGC 5291 as the disrupter without prejudging the issue of 
their origin). 

NGC 5291 itself may show the beginnings of a faint op- 
tical “tail” and dust lane extending to the north and slightly 
west, opposite to the position of the Seashell, altf ">ugh it is 
just as likely that these features may instead be part of a 
spiral arm. If they are tidal in origin however (i.e., it is a 
“tail” and not an “arm”), this type of morphology, delin- 
eated by tidal arms opposite a companion, with long patchy 
tails and scattered clumpy debris resembling dwarf galaxies, 
could be consist with a high velocity grazing interaction 
(Elmegreen et al 1993). 

Most of the features of the H l distribution detailed in 
these observations could probably be explained in terms of 
an interaction. Over three-quarters of the total H l mass is 
found in the extensions. The Seashell either is a gas free 
system or has been stripped of its atomic hydrogen, and there 
is more gas associated with the H f Concentration and with 
Knot B than with NGC 5291 itself. These characteristics are 
similar to those seen in NGC 7252 (Atoms-for-Peace) and 
NGC 3921, as reported by Hibbard & van Gorkom (1996). 
They examined a sequence of interacting/merging galaxies, 
and found that the more evolved the interaction, the lower 
the percentage of H I associated with the bodies of the parent 
galaxies. They also reported substantial offsets between the 
H l and the optical structures of their systems, and found that 
in the latest stage, essentially all the atomic gas w-as located 
in the tidal extensions. 

In addition to the optical and H l morphologies, the one- 
third solar metallicity of the gas in the knots is consistent 
w ith that expected in material stripped from the outer disk of 
a gas-rich system with an active star-forming history (al- 
though note the caveat about self-enrichment mentioned 
above). The H ll regions in the tails of the Antennae (NGC 
4038/39) also have metallicities about one-third solar (Mira- 
bel etal 1994). 

The kinematics of the gas (Figs. 7 and 8) may also be 
consistent with an interaction-dominated system. Interactions 
commonly produce plumes with monotonic slowly- 
increasing velocities away from the main bodies, especially 
if the material in the plumes is escaping, so that the outer- 
most features are not returning to the systemic velocity as 
occurs in 4 ‘bound” plumes like those in NGC 7252 (Fig. 4 
of Hibbard et al 1994). In an interaction, orbital decay oc- 
curs as energy is transferred from the galaxies to their par- 
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tides. The particles receive energy and angular momentum 
from the relative motion of the passage, which produces the 
tail-like structures and orbital decay. The tails produced by 
these encounters move radially away and fan outwards from 
the parent systems and the velocity increases along the tail. 
At the base, the tail particles are more tightly bound, turn 
quickly in their orbits and fall back towards the remnants. 
The more distant matter is less tightly bound and falls back 
more slowly, if at all. This monotonic increase in velocity 
along the tail is a prominent feature of the H I structure in the 
NGC 5291 system. The relative velocity reaches a maxi- 
mum, and then as the feature arcs to the NW and SW, the 
line-of-sight velocity component decreases as the tangential 
motions become more important. The velocity map (Fig. 8), 
for example, quite closely resembles that of the merger NGC 
4676 in Fig. 13 of Hibbard & van Gorkom (1996). In NGC 
5291 the reducing velocity differences of the material in the 
western arc could also perhaps be ascribed to the slowing of 
bound material as in NGC 7252, although it is also consistent 
with the velocity contours which might be expected if the 
western arc were an H I feature delineating the edge of a 
distorted disk. 

We can draw some additional conclusions about the is- 
sues to be considered if the interaction scenario is to be veri- 
fied. As noted, from the optical image the most obvious can- 
didate for the partner in this destructive minuet, owing to its 
contorted shape, is the Seashell. However, as first noted by 
L79, it is difficult to support this model if the passage is 
retrograde. The construction of tidal tails is most pronounced 
if the victims rotate in the same direction as the passage, and 
is further enhanced if the latter is roughly coplanar with the 
disk of the victim (i.e., the closer the alignment of the same- 
sign angular momentum vectors of the participants, the 
greater their disruption of one another). Such a prograde (but 
not necessarily coplanar) passage is needed to produce the 
morphology we observe in the NGC 5291 system. However, 
the Seashell appears to be approaching the observer in the 
region for which the H I in the system is receding. Therefore, 
a prograde orbit requires a special orientation for the NGC 
5291 -Seashell interaction, and one with a high velocity pas- 
sage. It is clear that careful modeling (along the lines of that 
discussed by Barnes 1994 and Hibbard & Mihos 1995) will 
be required to clarify the possibility of constructing a suc- 
cessful scenario. 

Alternatively, the Seashell may not be the culprit here. In 
addition to the orthogonal interaction/ring-galaxy hypothesis 
put forth by L79, M. Balcells (1996, personal communica- 
tion) has suggested that Knot B is an irregular that has en- 
tered the system at almost the escape velocity and is now 
spiraling into the center of NGC 5291; perhaps it was a sat- 
ellite to the Seashell. Such a model has many resemblances 
to the relation between M31, the Milky Way, and the LMC 
(Lynden-Bell 1994). The dense Hi distribution shown in 
Fig. 3 looks very much like a trail of debris torn from a 
system dominated by the Knot B complex. Much of the ring- 
arc might denote the orbit of the passage, although the 
amount of H I involved is considerable, possibly one-half the 
total mass. The projected difference in velocity between 
Knot B and NGC 5291 is 185 km s“ ! and the projected 


separation is 61 kpc. This would be a bound system if the 
mass of NGC 5291 were several times its calculated (mini- 
mum) value (Sec. 4.2). However, even one deeply plunging 
passage could be extremely damaging if not totally destruc- 
tive to the “Magellanic” system. 

If, however, this system has been produced by an interac- 
tion between NGC 5291 and the Seashell, some of the spe- 
cifics can be inferred from the velocity field. The relative 
velocity between NGC 5291 and the Seashell is either ap- 
proximately 450 km s~ ! (Due 1995) or 630 km s -1 (L79), 
and the projected separation is 10.7 kpc (L79). Using the 
lower value, the mass necessary to gravitationally bind the 
Seashell and NGC 5291 is then 5 X 10 m A/q . The indicative 
(keplerian) mass for NGC 5291 is 1.5X10 ll Af© (see Sec. 
3.2 , below), which is a factor of 3 less than this binding 
mass. However, there are the usual uncertainties associated 
with the unknown orbital phase terms and inclinations to the 
line of sight which easily bracket this figure. Furthermore, 
Barnes (1994) shows that the overall dynamics of encounters 
are largely controlled by the extended dark halos of the sys- 
tems. As part of the interaction, dark halos transfer angular 
momentum and energy, imparting spin and raising broad 
tidal features. It is this halo interaction that brings about the 
merger and extracts luminous tails from galactic disks. In the 
models considered by Barnes, the ratio of halo to disk 
masses is 5:1, which is within the range of corrections to 
visible masses found by van Moorsel (1987) in a study of 
binary galaxies. Given all the uncertainties, we believe that it 
is entirely possible that the Seashell and NGC 5291 are 
bound and merging. 

While somewhat suggestive of the debris fields produced 
in some galaxy encounters, the optical image (Fig. 1) alerts 
us to the unusual nature of this system. Although appearance 
is dependent on both geometry and viewing angle, the com- 
plex shoals of objects (known to be H II regions) to the north 
and south of the main body in Fig. 1 do not closely resemble 
the tails and countertails seen in canonical interacting or 
merging galaxies such as the “Antennae” NGC 4038/9 and 
the “Atoms for Peace” galaxy NGC 7252 (e.g., Hibbard & 
van Gorkom 1996). Nor, indeed, does this system closely 
resemble any of the systems in the atlases of Arp (1966) or 
Arp & Madore (1987). 

The classical interactors produce tails and countertails 
which are relatively narrow and quite smooth, sometimes 
studded along their lengths “string-of- pearls” fashion with a 
number of condensations, often H II regions, which generally 
number a dozen or less and are quite inconspicuous com- 
pared to the continuous structure. In the NGC 5291 complex, 
in contrast, there is little continuous structure discemable 
beyond the visible galaxy. There are perhaps two faint fea- 
tures extending ~1'5 to the SSE and, reinforced by an inner 
dust lane to resemble a weak spiral arm, slightly west of 
north (Fig. 1), but these do not underlie the main swarms of 
H II knots. It may also be noted that the H II knots in NGC 
5291 are mostly not arranged in a line (except, perhaps, to 
the south between D and G), but appear to be spread out 
laterally (i.e., knots H through M: the western arc) much 
more than is usual in the long filamentary arms and tails in 
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the classical interactors. A final discrepancy is evident in the 
unusually large number of the H n knots. 

The energetics of the gas in the NGC 5291 system imply 
that a significant fraction of the gas is unbound (at 20 kpc 
(72")) unless the galaxy contains a massive halo. The escape 
velocity is about 255 km s~ l , well below the maximum ve- 
locities observed). However, large velocities and energetics 
are apparently not rare in the early phases of galaxy mergers. 
Similar situations exist for NGC 4676 and especially for Arp 
295 (Hibbard & van Gorkom 1996). For these three systems 
(NGC 5291, NGC 4676, and Arp 295) a mass of order 
10 12 M o is required to bind the H I. Thus, either substantial 
amounts of gas escape from the systems or they possess dark 
matter. We presume both to be true. Certainly Hibbard & 
van Gorkom (1996) note the role of escaping gas in the sce- 
nario they describe for galaxy mergers. The existence of 
large kinetic energies in NGC 5291 is thus neither unusual 
nor unexpected, if in fact it is in the early stages of a galaxy 
merger. Mirabel and colleagues (Mirabel et cil 1991; Mira- 
bel et ai 1992), find similar situations for the Superantennae 
(with a tail 260 kpc in length) and Arp 105; at the tip of a 
tidal tail, star formation is occurring. These are objects with 
masses in excess of the LMC, and which Mirabel and his 
associates believe were detached from the progenitor galax- 
ies ~10 8 years ago. Star formation is now occurring well 
after ejection of the tidal debris. Knot B represents an analo- 
gous situation; we will argue that Knot B is a stable entity 
(Sec. 3.2). 

The formation and evolution of dwarf irregular galaxies 
from tidally removed material in interacting systems has 
been postulated by Zwicky (1956) and Schweizer (1978). 
Much attention has recently been devoted to clumps in tidal 
tails, as prominent Ho\ H t and stellar enhancements have 
been shown to exist in these regions (e.g., Mirabel et ai 
1991; Mirabel et al. 1992; Hibbard et al. 1994; Hibbard & 
van Gorkom 1996). This discovery of regions of star forma- 
tion associated with the tidal tails of interacting galaxies is 
evidence that gravitational condensation is possible in tidally 
removed material. Previous high spatial resolution observa- 
tions of the H I in tidal tails have shown that the distribution 
of gas in these regions is often clumped into organized re- 
gions of high density with masses on the order of lO g M 0 
(Hibbard et al. 1994). These concentrations are very similar 
to dwarf irregular galaxies in terms of their sizes, masses, 
colors, and stellar and gas distributions. The observed star 
formation rates are also very similar to those found in ac- 
tively star-forming dwarf irregular galaxies. 

Numerical modeling simulations have also been per- 
formed to test the possibility that material ejected during a 
tidal encounter will remain bound, or develop into sub-units 
on the scale of dwarf- and even full-sized irregular galaxies 
(e.g., Barnes & Hernquist 1992; Elmegreen et ai 1993). The 
models by Barnes & Hernquist (1992) produced strings of 
bound subsystems along the unbound tidal tails of two model 
spirals passing each other in a prograde, direct, initially para- 
bolic encounter prior to merging. They find in their simula- 
tions that distinct clumps form in tidal tails, and that these 
objects are gravitationally bound and stable against tidal 
stripping. These clumps form from the stellar compo- 


nent; the most massive of these clumps are able to capture 
significant amounts of gas. 

Rather similar results were obtained by Elmegreen et al. 
(1993). They claim that the interaction will increase the en- 
ergy of the interstellar medium (ISM) and thereby raise the 
Jeans mass. This process erases fluctuations smaller than 
10 8 Mq . They show that the formation of these objects at the 
ends of tidal tails is a natural consequence of the model, and 
that they will continue to exist as distinct objects long after 
the tails have faded. 

Hibbard et ai (1994), amongst others, suggest that the 
giant H II regions seen near the ends of the tidal tails in NGC 
7252 are examples of such objects; other candidates are dis- 
cussed by Mirabel and collaborators (Mirabel et ai 1991; 
Mirabel et ai 1992; Mirabel et ai 1994). Indeed, Hibbard & 
Mihos (1995), in their model of NGC 7252, find that at least 
one of these stable clumps will ultimately orbit the remnant 
as a satellite, with peri/apocentric distances of 14 and 120 
kpc and an orbital period of about 4 Gyr. 

It seems clear that interactive processes can an.i do pro- 
duce bound systems resembling dwarf irregulars in this way. 
However, objects like Knot B, or the features in the Super- 
antennae and Arp 105 discussed by Mirabel et ai (1991, 
1992), have gas masses that are an order of magnitude 
greater than the total masses of the clumps that appear in 
models. 

Since ram-sweeping is the most obvious explanation for 
the asymmetry of the atomic gas and the large x-ray halo and 
tail, it is possible that we are seeing an interacting system 
that is subsequently being ram-swept as it passes through an 
intracluster medium. However, neither the ram-swept large 
disk nor the ram-swept interaction hypothesis alone can eas- 
ily explain all of the observed characteristics of this unusual 
system. Detailed dynamical modeling is essential if more 
certain diagnoses are to be obtained. 

4.2 Analysis of Dwarf Galaxy Candidates in the NGC 5291 
System 

Regardless of the formation mechanism, it is clear that 
there are large star-forming regions associated with the NGC 
5291 system. As pointed out by L79, the optical knots in the 
NGC 5291 complex resemble dwarf galaxies in size and con- 
tent. Two conditions must be met for a concentration of ma- 
terial to remain bound in the outer regions of such a complex 
field. These conditions require that the object (1) remains 
bound against its own internal kinetic energy (is self- 
gravitating), and (2) remains bound against the tidal force ot 
the parent galaxy. We here examine whether the knots ob- 
served in the NGC 5291 complex can remain bound against 
these two disruptive influences, as then the likelihood of 
their eventual detachment from NGC 5291 to become inde- 
pendent systems is quite high (see L79). We note that the 
following arguments are independent of the process which 
has given rise to the overall structure of the complex. (Ram- 
sweeping in fact causes compression, providing extra con- 
tainment and facilitating the production of a bound sub- 
system; if this is the initiating process our estimates, below, 
will be even more conservative.) 
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To begin to understand the NGC 5291 system it is neces- 
sary to develop an estimate of the total mass present. This 
may be done in a number of ways, making different assump- 
tions and therefore arriving at substantially different results. 
The simplest of these is to assume, with L79, that the system 
is a bound rotating disk extending out to the outermost H I 
and H II features, and that their velocities are projected or- 
bital velocities. This results in mass estimates similar to 
those of L79 (who made this assumption), scaled for our 
larger assumed Hubble constant: M tot ~ 10 12 Mq , where we 
have implicitly assumed that the galaxy is not far from 
edgeways-on. This value is large, but not without precedent 
in early-type spirals: for example, NGC 5084 (Gottesman & 
Hawarden 1986) has a mass ^ 1 .5 X 1 0 12 Mq , though much 
less of this is H I than appears to be the case in NGC 5291. 

A second approach is to ignore the H I and H II features in 
the outskirts of the complex, as under the interaction hypoth- 
esis of their origin they may well not be bound to NGC 5291 
at all. A lower limit to the mass of NGC 5291 may be de- 
termined using a different Keplerian approximation: 

rv 2 

( 5 ) 


and estimating the orbital velocity as the full width at 20% 
peak intensity of the velocity profile of the H I in the central 
optical image, corrected for inclination: 


v — 


2& v 20 
sin i 


( 6 ) 


From Table 3 we have At; 2 o = 464 km s -1 . For the 
extent of the gas we adopt the optical major axis extent, l'l 
(NED) 6 and a distance of 58 Mpc, resulting in a linear diam- 
eter of 18.2 kpc. For an inclination of 60° (L79) the total 
mass becomes 1.5X Allowing for a dark halo, we 

increase this value by a factor of 5, as in the models by 
Barnes (1994). 


4.2.1 Binding mass 


For an object to maintain gravitational stability, its gravi- 
tational potential energy must equal or exceed its internal 
kinetic energy. To estimate the binding mass, M b , necessary 
for an object to meet condition (1), we start with the well- 
known energy equation for two gravitating point masses: 


i> 2 = G(M + m) 


2 

r 


1 

a ’ 


(7) 


where a is the semi-major axis of the orbit, and m is the 
mass of a “test particle’’ at a distance r. Setting a equal to 
infinity for a parabolic orbit and assuming M>m , we end up 
with the classic escape velocity equation. Rearranging, we 
find an expression for the minimum mass required to bind 
the system: 

rv 2 

(«) 


NASA/IPAC Extragalactic Database. 


We approximate the internal kinetic energy by means of the 
half-width of the H I velocity profile at the 20 level ( jAl» 2 o)- 
Assuming an isotropic velocity dispersion: 

v=j3CiAv 20 ). (9) 

The velocity profiles for NGC 5291 and the individual knots 
are shown in Fig. 11. Alternatively, if the H I is in an in- 
clined plane, the >/3 correction is equivalent to an inclination 
of 35°. Our final expression for the binding mass now be- 
comes 


M b = 


r [ V 3 ( 2^ v 20)] 2 

2G 


( 10 ) 


As long as the knots are essentially spherical, and we are 
only considering the dense H I cores, this is not likely to be 
much in error. 

We have evaluated the binding masses for the major knots 
in the NGC 5291 complex. To do this, we have estimated the 
radius of each object from the cores of the H I associated 
with each knot as delineated by the high resolution CnB 
array data. Knots D, J, and L have cores which are unre- 
solved by the CnB array beam, while F-G is nearly so. We 
have therefore used the beamwidth as an upper limit to their 
diameters. The velocity widths and H I masses have been 
determined from the spectra in Fig. 1 1 . These spectra were 
calculated by summing the flux in each channel of the more 
sensitive (but lower resolution) DnC array data cube, within 
the high density regions (“cores”) as delineated by the high 
resolution CnB array data; the method makes the exact 
shapes of the spectra somewhat uncertain. The velocity 
widths were measured at 20% of the peak value to eliminate 
some of the uncertainty due to these fluctuations. Table 3 
lists the velocity widths, radii, resultant binding masses, and 
H I masses for each object. 

To determine whether the knots were self-gravitating, we 
have compared the measured H I mass for each knot (Table 
3) to the binding mass. Values range from 0.82 for Knot B to 
0.02 for Knot C. However, this is a conservative test, which 
assumes that there is no other contributor to the mass than 
the observed H I. In the case of knots B and F-G there are 
clearly two other components: the stars and the ionized com- 
ponent in the optically-visible SFR associated with the H I 
feature, and observed in the radio continuum. Allowing for 
these masses, Knot B is self-gravitating, as is region F-G. 
Knots H, K and the H I Concentration may be marginally 
bound as well. The case for the stability of Knot B against its 
internal energy is clearly a very strong one. We would have 
to be underestimating the mass by more than a factor of 4 
before the properties of Knot B are peculiar compared to a 
typical irregular galaxy. 


4.2.2 Tidal condition 

To evaluate the second condition, we calculated the tidal 
limit for each dwarf candidate. The tidal radius is the radius 
within which material in a knot must lie to remain bound 
against the tidal force between the object and the parent gal- 
axy. An object whose actual radius is smaller than its tidal 
radius will be stable against tidal disruption. 

The tidally limited radius, r T> can be expressed by (von 
Hoemer 1957): 
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Fig 1 1 (a) Velocity profiles for components of NGC 5291 complex. Vertical axis (flux in Janskysl is drawn to same scale for all six objects. The peak 
intensity is 0.025 Jy. (b) Velocity profiles for components of NGC 5291 complex. Vertical axis (flux in Janskys) is drawn to same scale for all six objects. The 
peak intensity is 0.0125 Jy; note that the vertical scale is one-half of that used for (a). 
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Fig. 1 1 . (continued) 

m | 1/3 lite galaxy, M is the mass of the parent galaxy, m is the mass 

T G 3M ' ' ' of the satellite galaxy (both assumed to be point sources), 

and R g is the projected separation between the two galaxies, 
where r r is the tidal limit of a test particle orbit in the satel- Determining the tidal limit requires that the total mass of 
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NGC 5291, the candidate’s separation from it, and the radius 
of each object be known. The effective tidal radius (r T ) has 
been calculated for each knot using its measured separation 
distance from NGC 5291 (listed in Table 3) and the knot's 
Hi mass. The total mass estimate used for NGC 5291 is 
IX 10 n A/ o * which includes the factor of 5 to allow for a 
dark halo. 

This calculation provides a conservative tidal limit which 
is correct in detail for a point source potential, and quite 
reasonable for this system. Furthermore, the mass of an iso- 
thermal, spherical halo increases directly as R Cj * Thus r r is 
proportional to (R G ) 2iy . As our value for M cannot be re- 
duced and R g is underestimated owing to projection effects, 
our assumed values for R G and M should lead to a realistic 
lower limit to the tidal radius r T . We assume that the knots, 
as they are small, are essentially spherically symmetric and 
any model dependence on M propagates slowly at the one- 
third power. Therefore, we believe that the tidal radius can 
only be larger than indicated by Eq. (1 1). 

This tidal limit is compared to each object’s radius as 
measured from the H I data to determine its tidal stability. As 
seen in Table 3, most of the high-density cores of the objects 
meet this criterion. Knots C, M, and the H 1 Concentration do 
not, although for the H I Concentration, the H I radius ( 1 .9 
kpc) just barely exceeds the tidal limit (1.8 kpc). 

Candidate components that meet both the binding and 
tidal criteria are very possibly irregular galaxies in various 
stages of evolution. It may be argued that a subset of these 
objects become independent irregulars. An understanding of 
the likelihood of this process of galaxy formation will allow' 
us to estimate its contribution to the overall dwarf popula- 
tion, especially in clusters where interactions and ram- 
stripping are more likely to occur. 

5. CONCLUSIONS AND SUMMARY 

The high resolution H I observations presented and dis- 
cussed above shed considerable light on the nature of the 
NGC 5291 system, but do not finally resolve the causes of its 
peculiarities. Nevertheless we now' have a much clearer pic- 
ture of the basic parameters. 

NGC 5291 is being significantly perturbed by interaction, 
either with another galaxy, or with the intracluster medium 
of Abell 3574, or both. What remains unclear is the relative 
weight of these processes in arriving at the present condition 
of the system. However, we are at least in the position of 
being able to describe that condition in some detail. NGC 
5291 is attended by veritable shoals of star-forming H II 
complexes, the brighter specimens of which are associated 
with concentrations of H I on our VLA images. 

The H 1 observations establish clearly that a major specu- 
lation of L79 was in several elements correct: at least one 
(Knot B) and possibly several (H, K, and the H I Concentra- 
tion) of the H i/H II features are bound against their own 
internal energies and stable against tidal disruption by the 
nearby galaxy. These objects are therefore prime candidates 
for young irregular galaxies produced from stripped galactic 
material. An argument particularly may be made for Knot B, 
which possesses a large young stellar component (Due 1995) 


in addition to its H I mass and performance on the stability 
criteria. Knot B appears in every way to resemble an LMC- 
sized (or larger) “dwarf” irregular galaxy and would cer- 
tainly be categorized as such were it to be detached from its 
parent galaxy, as appears likely in the next gigayear or so: it 
is a true, young, galaxy. 

We suggest that the other knots may be dwarf irregular 
galaxies in formative stages of evolution; Knot B represent- 
ing the largest and perhaps most highly evolved of these 
candidates. An important implication of this work is the veri- 
fication that genuinely young galaxies may evolve from 
components of the debris fields of destructive interactions 
between parent systems and both their companions (via tidal 
effects) and their environments (via ram-sweeping). 

Understanding of the processes involved in this forma- 
tion, and hence, for example, the age of the new entity, re- 
quires solving the puzzle of the dominant process in the 
NGC 5291 system. We have attempted to set out the obvious 
evidence for the ram-sweeping of a large disk or of a tidally 
disrupted system, without selecting either. This is because 
we consider it essential that careful modeling of the system 
and its processes on both these scenarios be undertaken, if a 
correct understanding is to be obtained. Such an understand- 
ing is in turn essential if the lessons to be learned from the 
availability of a young dwarf irregular galaxy are properly to 
be absorbed. 

It has been suggested that halt or more of the dwarf ir- 
regular galaxies in compact clusters have formed from gravi- 
tational encounters of the more massive galaxies in the clus- 
ter (Hunsberger et al 1996). The observations presented 
herein, along with a growing body of evidence presented by 
many other researchers, seriously weaken models of galactic 
evolution that attempt to explain the various types of galax- 
ies seen in the universe as the result of different, independent 
processes. To date however, too few systems have been ob- 
served to determine if a new class of genuinely young gal- 
axies is truly indicated, or if instead this formation process is 
rare in the universe. Further investigations are needed to es- 
tablish the extent to which gravitational encounters are in- 
volved in the building of dwarf irregular galaxies. 
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lex taken with the ESO 3.6 m telescope (Pedersen et al. 1978: reproduced by permission), 
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